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Abstract 

A search for squarks in i?-parity violating supersymmetry is performed in e^p colh- 
sions at HERA using the HI detector. The data were taken at a centre-of-mass energy of 
319 GeV and correspond to an integrated luminosity of 64.3 pb~^ for e~^p coHisions and 
13.5 pb~^ for e~p collisions. The resonant production of squarks via a Yukawa coupling 
A' is considered, taking into account direct and indirect i?-parity violating decay modes. 
No evidence for squark production is found in the multi-lepton and multi-jet final state 
topologies investigated. Mass dependent limits on A' are obtained in the framework of the 
Minimal Supersymmetric Standard Model. In addition, the results are interpreted in terms 
of constraints on the parameters of the minimal Supergravity model. At the 95 % confi- 
dence level squarks of all flavours with masses up to 275 GeV are excluded in a large part 
of the parameter space for a Yukawa coupling of electromagnetic strength. For a coupling 
strength 100 times smaller, masses up to 220 GeV can be ruled out. 
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1 Introduction 



The ep collider HERA is ideally suited to search for new particles coupling to electron^ -quark 
pairs. In supersymmetric models (SUSY) with i?-parity violation (J/lp), squarks can couple to 
electrons and quarks via Yukawa couplings A'. At HERA, squarks could be produced resonantly 
via the fusion of the incoming 27.6 GeV electron and a quark from the incoming 920 GeV 
proton. Squark masses up to the electron-proton centre-of-mass energy, ^/s — 319 GeV, are 
kinematically accessible. 

This paper describes a search for squarks of all flavours using HI data corresponding to 
an integrated luminosity of 64.3 pb~^ for e+p collisions and 13.5 pb~^ for e~p collisions. The 
search is carried out in the framework of the Minimal Supersymmetric Standard Model (MSSM) 
in the presence of a non-vanishing A'. The analysis covers the major event topologies from di- 
rect and indirect squark decays such that the results can be interpreted in terms of a wide 
range of SUSY parameters. The search presented here supersedes the results previously ob- 
tained by HI [1, 2] at a lower centre-of-mass energy (^/s ?a 300 GeV) and with fewer data. 
Complementary searches for SUSY have been carried out at the LEP e+e~ collider [3, 4] 
and at the TeVatron pp collider [5, 6]. 



2 Phenomenology 

In the most general supersymmetric theory that is renormalisable and gauge invariant with re- 
spect to the Standard Model (SM) gauge group, the i?-parity Rp = (_i)3b+l+25^ where B 
denotes the baryon number, L the lepton number and S the spin of a particle, is not conserved. 
Couplings between two ordinary fermions and a squark (g) or a slepton (/) are then allowed. 
The Yukawa couplings responsible for squark production at HERA are described in the su- 
perpotential by the terms \'--j^LiQjDk, where i, j and k are family indices. Lj, Qj and are 
superfields, which contain the left-handed leptons, the left-handed quarks and the right-handed 
down quark, respectively, together with their SUSY partners and o?^. The corresponding 
part of the Lagrangian expanded in fields is given by 



^LiQjDu - Kjk 



+ C.C., (1) 



where the superscript c denotes the charge conjugate of a spinor and * the complex conjugate of 
a scalar field. Non-vanishing couplings A'^^;. allow the resonant production of squarks at HERA 
through eg fusion [7]. The values of the couplings are not fixed by the theory. For simplicity, it 
is assumed here that one of the A'^^;. dominates over all other possible trilinear couplings. 

For the nine possible couplings A'^^^;., the corresponding squark production processes in e^p 
reactions are listed in table 1. At high Bjorken-x the density of antiquarks in the proton is 
smaller than that of quarks. Thus e^p scattering gives the best sensitivity to the couplings 

'in the following the term electron will be used to refer to both electron and positron unless explicitly stated 
otherwise. 
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Table 1: Squark production processes in e^p collisions for different Yukawa couplings A'^^^. 
The qr and ql symbols denote the squarks which are superpartners of the right- and left-handed 
quarks, respectively. Their antimatter counterparts are denoted by qj^ and g^. 



'^'iifc = 1)2,3), where mainly d^-type squarks (i.e. the superpartners dR,SR and bR of 
right-handed quarks) can be produced. The dominant squark production cross section in e~p 
collisions is approximately proportional to X'^^;^ ■ u{x) where u{x) is the probability of finding 
a u quark in the proton with a momentum fraction x = M?/s, Mq being the squark mass. 
In contrast, e^p scattering gives the best sensitivity to the couplings A'^j^ (j = 1, 2, 3), where 
mainly CtL-type squarks {i.e. the superpartners ul,cl and of right-handed quarks) can be 
produced. Here the dominant squark production cross section is approximately proportional to 
A']^^ ■ d{x). Since the u quark density in the proton is larger than the d quark density at large 
,T, the squark production cross section in e~p interactions is larger than that in the e^p case for 
comparable couplings . 

In this work the signal cross section is obtained in the narrow width approximation from the 
leading order (LO) amplitudes given in [8], corrected by multiplicative factors [9] to account 
for next-to-leading order QCD corrections. The parton densities are evaluated at the hard scale 
M|. For cases, where the squark width is not negligible, the approach given in [10] is followed. 

In SUSY all supersymmetric particles are unstable. Squarks can decay via their Yukawa 




Figure 1 : Lowest order s-channel diagrams for squark production via the Yukawa coupling 
A' in (a) e~p and (b) e+p interactions, followed by squark decays. 
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coupling A' into SM fermions. According to equation (1), the d/j-type squarks can decay either 
into e~ + vP or z/g + d\ while the tt^^-type squarks decay into e+ + only. The squark 
decays, proceeding directly via the couplings A'^^^^ and A'^^-^, are illustrated in figure 1. 

Squarks can also decay via their usual Rp conserving gauge couplings, as shown in figure 2. 
The Mz,-type squarks can undergo a gauge decay into states^ involving a neutralino x° (i ~ 
1 ... 4), a chargino (i = 1, 2) or a gluino g. In contrast, d^-type squarks decay to Xi or 
g only and decays into charginos are suppressed, since the supersymmetric partners of right- 
handed quarks do not couple to winos. 

The final state of these gauge decays depends on the subsequent gaugino decay, of which 
examples are shown in figure 2. Neutralinos Xi with i > 1 as well as charginos (gluinos) are 
expected to undergo gauge decays into a lighter x and two SM fermions (two quarks), through 
a real or virtual gauge boson or sfermion (squark). The decay chain ends with the ^p decay of 
one sparticle, usually the lightest supersymmetric particle (LSP), assumed here to be a x'^, x"^ 
or g. l^p decays of gauginos are mainly relevant for the lightest states. Neutralinos may undergo 
the ]^p decays x° e^qq' or x° vqq, the former (latter) being more frequent if the x° is 
dominated by its photino (zino) component. Gluinos can undergo the same I^p decays. When 
a x° or a ^ decays via an ]^p coupling into a charged lepton, both the "right" and the "wrong" 
charge lepton (with respect to the incident beam) are equally probable, the latter case leading to 
striking, largely background-free signatures for lepton number violation. In contrast, the only 
possible decays for charginos are x"^ — vu^d^ and x^ ~^ e^d^d^. 

The u'j^ (dji) decay chains analysed in this paper are classified by event topology, as de- 
scribed in table 2. This classification relies on the number of charged leptons and/or hadronic 
jets in the final state, and on the presence of missing momentum. The channels labelled eq and 
vq are the squark decay modes which proceed directly via ^p couplings, while the remaining 
channels result from the gauge decays of the squark and are characterised by multijet (MJ) final 
states. The channels labelled e+M7, e~MJ and uMJ involve one or two SUSY fermions (x 
or g) denoted by X and Y in table 2. The channels e£MJ and uiMJ necessarily involve two 
SUSY fermions. 

^The mass eigenstates Xi (xf) ^re mixed states of the photino, the zino and the neutral higgsinos (the winos 
and the charged higgsinos). The g is the SUSY partner of the gluon. 
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Decay patterns involving more than two gauginos are kinematically suppressed and are 

therefore not expHcitly studied here. Processes leading to final states with tau leptons are also 
not explicitly investigated. Cases where the Xi has such a long lifetime that large displaced 
vertices are expected are not considered, since the region of parameter space that allows a Xi 
to escape detection for a finite value of the ]/lp coupling is very strongly constrained by the 
searches for gauginos carried out at LEP [4]. Decays of a x into states involving a Higgs boson 
are taken into account when the Higgs decays into hadrons. The contribution of these decays 
is, however, very small. 



3 The HI Detector 

A detailed description of the HI experiment can be found in [1 1]. The main components of the 
tracking system are the central drift and proportional chambers which cover the polar angle^ 
range 25° < 9 < 155°, a forward track detector (7° < 9 < 25°) and a backward drift chamber. 
The tracking system is surrounded by a finely segmented liquid argon (LAr) calorimeter [12] 
which covers the polar angle range 4° < < 154° and has an energy resolution of a{E)/E ~ 
12%/^E/GeV © 1% for electrons and a{E)/E ~ 50%/^/E/GeV © 2% for hadrons, as 
obtained in test beam measurements [13]. The tracking chambers and the LAr are surrounded 
by a superconducting solenoid and its iron yoke instrumented with streamer tubes. The latter 
are used to detect hadronic showers which extend beyond the LAr and to identify muons. The 
luminosity is determined from the rate of Bethe-Heitler events (ep — > ep7) measured in a 
luminosity monitor. 



4 Monte Carlo Event Generation 

In order to estimate the amount of SM background in the various squark decay channels and 
to determine the signal detection efficiencies, complete simulations of the HI detector response 
are performed for various Monte Carlo (MC) samples. 

For each possible SM background source, a sample of MC events is used, corresponding 
to a luminosity of more than ten times that of the data. The determination of the contribution 
of neutral current (NC) deep inelastic scattering (DIS) processes is performed using two MC 
programs which both include LO QCD matrix elements but employ different models of QCD 
radiation. The first is produced with the DJANGO [14] event generator, where QCD radiation 
is implemented using ARLADNE [15], based on the Colour Dipole Model (CDM) [16]. This 
sample is chosen to estimate the NC DIS contribution in the eq channel. The second sample 
is generated with the program RAPGAP [17], where higher order QCD radiation is modelled 
using leading-log DGLAP parton showers [18]. This sample is used to determine the NC DIS 
background in the final states with an electron and multiple jets, as RAPGAP gives the better 
description of this particular phase space domain [10]. For both samples, the parton densities 
in the proton are taken from the CTEQ5L [19] parameterisation. Hadronisation is performed 

^The polar angle 6 is measured with respect to the direction of the outgoing proton beam i+z). 
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Channel 


Decay process 


Event topology 


eq 


A' 

q — > e q 


high pt e + I jet 


uq 


a|j — > Ve a 


missing px + 1 jet 


e^MJ 


q ^ q X 

^ e^qq 
q ^ q X 

^ qq Y 

^ e^qq 


e (both charges) 
+ multiple jets 


uMJ 


q ^ q X 

A' 

lyqq 

q — > q X 

qq Y 
A' 

^ lyqq 

q ^ q X 

^ up Y 

A' 

^ vqq 


missing pt 
+ multiple jets 


elMJ 


q ^ q X 

iue Y 

^' ±- 
^ e qq 

q — > q X 

^ e+e- Y 

^' ±- 

^ e qq 

q ^ q X 

^ e+e- Y 

A' 

^ uqq 


e 

+ £ (e or yu) 
+ multiple jets 


ulMJ 


q ^ q X 

eue Y 

A' 

^ uqq 

q > q X 

^ up Y 

A' 

^ eqq 

q ^ q X 

li^li- Y 

A' 


i (e or /i) 
+ missing pt 
+ multiple jets 



Table 2: Squark decay channels in SUSY classified by event topology. X and Y denote a 
neutralino, a chargino or a gluino. The process is indicated by A'. 
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in the Lund string fragmentation scheme using JETSET [20]. The modelling of the charged 

current (CC) DIS process is performed using the DJANGO program with CTEQ5L parton den- 
sities. The direct and resolved photoproduction (7^) of light and heavy flavours, including 
prompt photon production, is generated using the PYTHIA [21] program, which relies on first 
order matrix elements and uses leading-log parton showers and string fragmentation. The SM 
expectations for ep — > eW^X and ep — > eZ^X are calculated using EPVEC [22]. The LO 
MC simulations used to model QCD multi-jet production give only approximate descriptions 
of the kinematic distributions. From the comparison of the distributions of multi-jet events be- 
tween the data and the LO MC simulations a normalisation factor of 1.2 is derived [10] which 
is applied to the yield of multi-jet events predicted by RAPGAP and PYTHIA. 

To allow a model independent interpretation of the results, all squark decay processes given 
in table 2 are simulated separately for a wide range of masses of the SUSY particles involved. 
The LEGO [23] event generator is used for the determination of the signal detection efficien- 
cies in the eq and uq channels, whereas for the gauge decays of squarks the SUSYGEN [24] 
generator is used. The squark mass is varied from 100 GcV to 290 GeV in steps of typically 
25 GeV. For gauge decays of squarks involving a gaugino which decays directly via (i.e. 
processes corresponding to the first line of the e^MJ and vMJ rows in table 2), the process 
q — > g'x5 is simulated for x5 masses ranging between 30 GeV and Mq. In order to study the 
cascade gauge decays which involve two gauginos, the processes q — > qxt QXiff' ^nd 
Q ~^ 1X2 ~^ QXiff' simulated for xt and X2 masses ranging between 40 GcV and Mg, and 
for Xi masses between 30 GeV and or M^o. The masses of the x's are varied in steps of 
typically 10 GeV. The lower mass values for squarks and x's are motivated by the exclusion 
domains resulting from SUSY searches at LEP [3,4]. The simulations allow the determina- 
tion of signal detection efficiencies as a function of the masses of the SUSY particles involved, 
since the mass intervals are sufficiently small for linear interpolations to be used. 



5 Searches for SUSY Signals 

5.1 Basic event selection 

The recording of the events used in this analysis is triggered using the LAr system [12], with 
an efficiency close to 100 %. Background events not related to ep collisions are suppressed by 
requiring that a primary interaction vertex be reconstructed within ±35 cm in z of the nominal 
vertex position and by using topological filters against cosmic and proton-beam related back- 
ground. The event time as determined by the central drift chambers is required to be consistent 
with the bunch crossing time. 

5.2 Particle identification and kinematic reconstruction 

The following criteria are used to select events containing leptons, high transverse momen- 
tum jets or missing transverse energy. An electron is identified as an isolated and compact 
electromagnetic cluster of energy greater than 11 GeV in the LAr. For electrons in the central 
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detector region (30° < 9^ < 145°) a charged track pointing to the electromagnetic cluster is 
required. A muon candidate is identified as a track measured in the central or forward tracking 
system, which matches geometrically with a track in the instrumented iron, a track in the for- 
ward muon detector or an energy deposit in the LAr calorimeter compatible with that expected 
from a minimum ionising particle. Hadronic jets are reconstructed from energy deposits in the 
calorimeter using a cone algorithm in the laboratory frame with a radius ^/Arf~+^^ — 1, 
where 77 = — In tan | is the pseudorapidity and (p denotes the azimuthal angle. The missing 
transverse momentum pT.miss is obtained by the summation of all the energy deposits in the 
calorimeter. 

For further selection the following Lorentz invariants are important: 

They are determined using the measurement of the polar angle Og, the energy and the trans- 
verse momentum p^.e of the electron with the highest pr found in the event. denotes the 
energy of the incident electron. Similar quantities can be calculated using the Jacquet-Blondel 
method [25]: 

where pT,h and Ylii^ ~Pz)h are calculated from the hadronic energy deposits in the calorimeter. 



5.3 Systematic uncertainties 

In each selection channel the systematic errors on the SM background expectation are evaluated 
by considering the following uncertainties. 

• The uncertainty on the electromagnetic energy scale of the calorimeter varies from 0.7 % 
to 3 % depending on the calorimeter region [26]. 

• The uncertainty on the hadronic energy scale is 2 %. 

• The uncertainty on the integrated luminosity is 1.5 %. 

• An uncertainty of ±7 % on the DIS expectation arises from the parton densities of the 
proton at high x. 

• An uncertainty of ±10 % on the predicted cross section for multi-jet final states is esti- 
mated by comparing the LO MC simulations where higher order QCD radiation is mod- 
elled by either the CDM or DGLAP parton showers. 

Furthermore, the following uncertainties related to the modelling of the SUSY signal are taken 
into account. 
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• The theoretical uncertainty on the signal cross sections due to uncertainties in the parton 
densities varies from 7 % for e~u d\ at low squark masses up to 50 % for e+d — > 

at high masses. 

• Choosing either or the square of the transverse momentum of the final state lepton in 
squark decays (proceeding directly via the coupling A') instead of Af| as the hard scale at 
which the parton distributions are determined yields an uncertainty of ±7 % on the signal 
cross section. 

• An uncertainty of 10 % is attributed to the signal detection efficiencies, resulting mainly 
from the interpolation between the simulated mass values. 

5.4 R-parity violating squark decays 
5.4.1 Selection channel eq 

The final state with an electron and a jet of high transverse momentum, resulting from squarks 
decaying in the channel eq, corresponds exactly to the NC DIS signature at high x. However, 
the Me and Ue distributions of the two processes differ. Squark decays via lead to a res- 
onance in the Mg distribution which is measured with a resolution of between 3 and 6 GeV 
depending on the squark mass. Squarks produced in the s-channel decay isotropically in their 
rest frame, leading to a flat da/dy distribution. In contrast, the distribution for NC DIS varies 
approximately as da/dy oc 

The selection criteria for the eq channel are the following. 

• The total transverse momentum of the events must be balanced: pr,miss < 15 GeV. 

• The reconstructed momentum loss in the direction of the momentum of the incoming 
electron must be such that AO < '^{E — p^) < 70 GeV, where the sum extends over all 
reconstructed particles. 

• An electron must be found in the LAr calorimeter withpT,e > 16 GeV. 

• To improve the sensitivity, the differences in the Me and ye distributions of the SUSY sig- 
nal and the DIS background are exploited by applying a lower yg-cut which depends on 
the mass of the squark under consideration. The ^e-cut is optimised by minimising the ex- 
pected limit. It ranges from 0.5 for masses around 100 GeV to 0.2 around 290 GeV [10]. 

• The selection is restricted to the kinematic range Ql > 2500 GeV^ and y^ < 0.9. Ex- 
cluding the highest y^. values avoids the region where migration effects due to initial state 
QED radiation are largest. Furthermore, background from photoproduction events, in 
which hadrons are misidentified as electrons, is suppressed. 

• To ensure that the various selections are exclusive"^, all events accepted in one of the 
selection channels with an electron and several jets (section 5.5) are not accepted in the 
eq channel. About 10 % of the candidate events are removed from the eq channel by this 
requirement. 
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Figure 3: Mass spectra for the eq selection channel in (a) e+p and (b) e^p collisions. The shaded 
error band indicates the uncertainty on the SM background. The signal expected for a squark 
of mass 150 GeV is shown with arbitrary normalisation (dashed histogram). Events with an 
electron and multiple jets are not included in the spectra. 

The Me spectra for the data and the SM background simulation after this selection are shown 
in figure 3 for e~^p and e^p collisions. No significant deviation from the SM expectation is 
found for either data sample. Table 3 gives the total numbers of selected events and the SM 
expectation. In the e^p data set 632 candidate events are found, which is to be compared with 
628 ± 46 expected from SM processes. In the e~p data sample, 204 events are observed while 
the SM expectation is 192 ± 14. 

5.4.2 Selection channel uq 

Squarks undergoing a direct |?p decay into uq lead to CC DIS-like events with high missing 
transverse momentum. The events are expected to cluster in the Mh distribution with a resolu- 
tion of 10 to 20 GeV, depending on the squark mass. 

The selection criteria for the uq channel are the following: 

• The missing transverse momentum must be greater than 30 GeV. 

• No electron or muon must be found with p-r > 5 GeV. 

• The events must lie in the kinematic range Q\ > 2500 GeV^ and i/h < 0.9. The resolu- 
tions in both Mh and Q\ degrades with increasing yh since both 6Mh/Mh and 5Q\/Q\ 
behave as 1/(1 — y/i) for yh ^ 1. Hence the high yh range is excluded. 

• To ensure exclusivity with respect to the uMJ channel (section 5.6), events with two or 
more jets with prjet > 15 GeV are rejected. This removes about 3.5 % of the candidate 
events. 

"^This is necessary for the limit calculations (section 6). 
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Figure 4: Mass spectrum for the uq selection channel in e^p collisions. The shaded error band 
indicates the uncertainty on the SM background. The signal expected for a squark of mass 
150 GeV is shown with arbitrary normalisation (dashed histogram). Events accepted in the 
selection channel vM J are not included in the spectra. 

Only rf/j-type squarks, which are produced mainly in e~p collisions, can undergo a decay into 
the uq final state. The Mh spectrum of this data set and the SM background are shown in 
figure 4. No significant deviation from the SM expectation is found. 261 events are observed in 
the data and 269 ± 21 are expected according to the SM. 

5.5 Squark gauge decays leading to e + jets + X final states 

For the channels e^MJ, e~MJ, eeMJ, efxMJ and ueMJ a common preselection is carried 
out: 

• At least one electron must be found with pT,e > 6 GeV in the angular range 5° < 6'e < 
110°. For central electrons {9e > 30°) the charged track, measured in the central tracking 
system, must geometrically and kinematically match the electromagnetic cluster. To dis- 
criminate against fake-electron background from photoproduction, electron candidates in 
the forward region {Of. < 30°) have to fulfill harsher isolation criteria and the ~ Pz) 
of the event must be greater than 30 GeV. The latter cut causes only a small efficiency 
loss for all channels discussed here. 

• At least two jets must be found with prjct > 15 GeV in the range 7° < ^^jct < 145°. 

• For all final state topologies considered here, the squark decay products are mainly emit- 
ted in the forward direction. This is exploited by requiring that: 

- Ql> 1000 GeV^. 

- At least one of the polar angles of the highest px electron and the two highest p^ 
jets is less than 40°. 

- Of the two jets with highest pt, that with the larger polar angle 6'backw satisfies 
6'backw < 180° ■ {ue — 0.3). This cut efficiently separates the SUSY signal events 
from the NC DIS background [10]. 
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Figure 5: Mass spectra for (a) the e+Af J selection channel in e+p collisions and (b) the MJ 
selection channel in e^p collisions. The shaded error band indicates the uncertainty on the 
SM background. The signal expected for a squark of mass 150 GeV is shown with arbitrary 
normalisation (dashed histogram). 

After this preselection, 91 (22) events are found in the e+p (e~p) data sample while 89.3 ± 3.7 
(22.6 ± 0.7) is the SM expectation. Further cuts are applied for each sub-channel. 



5.5.1 Channels with "wrong" and "right" lepton charge 

For the channels e^MJ and e~MJ no neutrinos are involved in the final state. Therefore the 
missing momentum is restricted according to j9r miss < 15 GeV and 40 < J2i^~Pz) < 70 GeV. 
To ensure that the selection is exclusive with respect to the eeMJ and ejjiMJ channels, events 
with a second electron with pT,e > 5 GeV and 5° < 6'e < 110°, or a muon with p^^^ > 5 GeV 
and 10° <9f,< 110°, are rejected. 

Events are accepted in the channel having the "wrong" charge lepton, i.e. different from 
the incident beam, if the electron/positron is found in the angular range 9e > 30° (where the 
charge measurement is made with the central tracking system) and the charge is measured to be 
opposite to that of the incident lepton, with a significance greater than two standard deviations. 
No candidates are found in the data and the SM expectation in this channel is very low (see 
table 3). 

In the "right" charge lepton channel, i.e. the same charge as the lepton beam, events are 
accepted if they contain either a central electron {9e > 30°) with a charge measurement of the 
"right" sign or an electron found in the forward region (9e < 30°). In the latter case no charge 
requirement is made. For the selected events, an invariant squark mass Afinv is calculated as: 
Minv = ^/4:E^ ~ -^e)' where the sum runs over the electrons and the jets found in the 

event with > 5 GeV. This method yields a good reconstruction of the squark mass with a 
typical resolution of 7 to 10 GeV. The M\^^ distributions for the data and the SM expectation are 
shown for the "right" charge eMJ channel in figure 5 for e+p and e~p collisions. No significant 
deviation from the SM is observed at any mass value. In total 72 (20) events are selected in the 
e+p (e^p) data set with 67.5 ± 9.5 (17.9 ± 2.4) predicted from SM background processes. 
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Figure 6: Mass spectra for the vMJ selection channel in (a) e+p and (b) e^p collisions. The 
shaded error band indicates the uncertainty on the SM background. The signal expected for a 
squark of mass 150 GeV is shown with arbitrary normalisation (dashed histogram). 

5.5.2 Channels with an additional lepton 

The further selection for the channels eeMJ and efxMJ requires either an additional electron 
with the same criteria as described in the common preselection, or an additional muon with 
Pt,^i > 5 GeV in the polar angle range 10° < 9^ < 110°. To ensure exclusivity, events are ac- 
cepted in one selection channel {eeMJ or efxMJ) only^. In these channels the SM background 
(mainly NC DIS) is very low (see table 3). No candidate events for either eeMJ or ejjiMJ are 
found in either data set, which is compatible with the SM expectation. 

For the channel veMJ, candidate events, possibly containing a neutrino, are selected by 
requiring pT,miss > 15 GeV. A cut of UeiVe — Uh) > 0.04 exploits the fact that for the SUSY 
signal the escaping neutrino carries a non-negligible part of J2 ~ Pz) and hence the variable 
Uh is substantially smaller than ye, while ye ~ yh is expected for NC DIS events. Events 
previously accepted in the eeMJ or ejjiMJ channels are rejected. No events are found in this 
channel. This is compatible with the SM expectation (mainly NC DIS) as detailed in table 3. 

5.6 Squark gauge decays leading iov+ jets + X final states 

The selection of vMJ and vfiMJ candidates starts with the requirement that: 

• The missing transverse momentum satisfies pT,miss > 26 GeV. 

• At least two jets with PTjet > 15 GeV are reconstructed in the angular range 7° < 6'jct < 
145°. 

• No electron is found in the event. 

''Events with a muon with px^^ > 5 GcV and 10° < 0^ < 110° are not accepted in the eeMJ channel. 
Similarly, events with an additional electron in the range 5° < 9e < 110° are not accepted in the efiMJ channel. 
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Table 3: Total numbers of selected events, SM expectations and ranges of selection efficiencies 
of the squark decay channels considered in e+p and in e~p collisions. The CtL-type squarks (e+p 
collisions) cannot decay to vq. 

If no muon is found, the event is identified as a vMJ candidate. Assuming that the missing 
energy of a candidate event is carried by one neutrino only, its kinematics are reconstructed 
by exploiting energy-momentum conservation. The four- vector of this v is then added to that 
of the hadronic final state to reconstruct the invariant mass Mrec of the incoming electron and 
quark. The mass resolution of this method is about 15 GeV. The Mrec spectra of the data and 
the expected SM background are shown in figure 6. In the e^p (e~p) data set 30 (12) uMJ 
candidate events are selected while 24.3 ± 3.6 (10.1 ± 1.4) are expected from SM background 
(mainly CC DIS). 

If in addition to the above requirements a muon with pt,^ > 5 GeV in the range 10° < < 
110° is found, the events are identified as vjiMJ candidates. No candidate events are found in 
either data set. This is compatible with the SM expectation (predominantly CC DIS), which is 
shown in table 3. 



6 Exclusion Limits 

The total numbers of selected and expected events are summarised in table 3 for all final state 
topologies considered in this analysis. It is assured by the choice of the selection cuts for the 
different channels that the selection of all topologies is fully exclusive. No significant devia- 
tion from the SM expectation is found in any channel. The selection channels are combined, 
separately for the e+p and e'p data sets, to derive constraints on SUSY models. 

6.1 Method of limit derivation 

For a given set of parameters in a certain supersymmetric model, the full supersymmetric mass 
spectrum and the branching ratios of all squark decay modes are calculated using the SUSY- 
GEN package. An upper limit iViim on the number of events coming from squark production is 
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calculated at a confidence level (CL) of 95 % using a modified frequentist approach based on 
Likelihood Ratios [27]. The following quantities enter the limit calculation. 

• The numbers of events observed in the data for all selection channels. For the channels 
where the SM background is considerable (eq, uq, eMJ with "right" lepton charge and 

uMJ) the event numbers are integrated within a mass bin around the squark mass under 
consideration. For each decay channel the width of the mass bin is optimised, i.e. the 
expected limit is minimised, using the reconstructed mass distributions from the SUSY 
signal and SM background simulations. For the channels eMJ with "wrong" lepton 
charge, eeMJ, ejiMJ, ueMJ and vjiMJ, no mass restriction is imposed, since the SM 
backgrounds are small. 

• The event numbers expected from SM background processes and their systematic uncer- 
tainties. 

• The signal detection efficiencies (see table 3) and their uncertainties for all squark decay 
processes in all selection channels, obtained using the calculated spectrum of sparticle 
masses. 

• The calculated branching ratios of all squark decay modes. 

A bound on the squark production cross section cxiim is then obtained from N\i^. Sets of model 
parameters that lead to signal cross sections above aiim can be excluded. 

The case of non- vanishing Yukawa couplings \'^^^ or A^^3, which correspond to the resonant 
production of stop and sbottom squarks, is treated separately since the top and bottom quark 
masses cannot be neglected in the calculation of couplings and branching ratios. Furthermore, 
a top quark could be produced in gauge decays. The top quark decays via f bW, leading to 
decay products different from those of the first two generations for which the efficiencies are 
determined. Diagrams which lead to a top in the final state are thus not taken into account in 
the calculation of the branching ratios. This represents a conservative approach, since most of 
the top decays are implicitly covered in the selection channels and would be visible in the mass 
distributions and the total event numbers. 

6.2 Limits in the "phenomenological" MSSM 

A version of the MSSM is considered here where the masses of the neutralinos, charginos and 
gluinos, as well as the couplings between any two SUSY particles and a SM fermion/boson, 
are determined by the usual parameters. These are the "mass" term which mixes the Higgs 
superfields, the SUSY soft-breaking mass parameters Mi, M2 and M3 for U{1), SU{2) and 
SU (3) gauginos, respectively, and the ratio tan (3 of the vacuum expectation values of the two 
neutral scalar Higgs fields. The parameters are defined at the electroweak scale. The gaugino 
mass terms are assumed to unify at a Grand Unification (GUT) scale to a common value mi/2, 
leading to the usual relations [28] between Mi, M2 and M3. The gluino mass is approximated 
by the value of M3 at the electroweak scale. The sfermion masses are free parameters in this 
model. Possible mixing between sfermions is neglected and all squarks are assumed to be 
degenerate in mass. The possibility of a photino-like Xi is first discussed, before turning to a 
complete scan of the SUSY parameter space. 
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Figure 7: (a, b) Exclusion limits at the 95 % CL on (a) A'^^^ with j = 1,2 and (b) A^^^;, with 
= 1, 2. (c - f) Branching ratios to the decay channels considered in the analysis for A' values 
at the exclusion limits shown in (a) and (b). The results are shown for MSSM parameters leading 
to a dominated by its photino component when slepton and squark masses are assumed to be 
degenerate (c, d) and for a slepton mass of 90 GeV (e, f). 
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6.2.1 Exclusion limits for a photino-like Xi 

For an example set of the MSSM parameters (// = -200 GeV, M2 = 80 GeV, tan/i = 2) 
leading to a Xi dominated by its photino component, exclusion limits at the 95 % CL on 
KjiU = 1' 2) and X[^k{k = 1, 2) are shown in figures 7 (a, b) as a function of the squark 
mass. The full curves represent cases in which sleptons and squarks are assumed to be degener- 
ate in mass. The dashed curves indicate the limits for slepton masses fixed at 90 GeV, close 
to the lowest mass bound from sfermion searches at LEP [3]. The HERA sensitivity allows 
tests of Yukawa couplings A' down to around 10^^ for squark masses of 100 GcV. For a high 
squark mass the sensitivity degrades since the production cross section decreases. At a squark 
mass of 290 GeV, A'^^^ (Kik) values larger than 0.6 (0.3) are ruled out. 

The branching ratios to all channels calculated for a A' value exactly at the exclusion limit 
are illustrated in figures 7 (c - f). The total branching fraction covered exceeds 75 % for all 
points in the MSSM parameter space and is generally close to 100 %. At large squark masses, 
a large Yukawa coupling A' is necessary to allow visible squark production. As a result the 
decay channels eq and uq proceeding directly via A' become important. For smaller masses, the 
dominant channels in the case of a photino-like Xi are those with an and several jets in the 
final state. 

For UL-type squarks (X'lji 7^ 0) the relative contributions of the gauge decay channels 
strongly depend on the slepton mass. In the case of a light slepton (Mj- — 90 GeV), the de- 
cays of a Xi iiito a lepton-slepton pair are kinematically allowed. Thus cascade gauge decays 
of M^-type squarks are possible, leading to enhanced contributions from the channels elMJ and 
ulMJ. In contrast, the cascade gauge decays of w^-type squarks are kinematically suppressed 
for Ml = Mg. The dependence of the X[j-^ limit on the slepton mass is rather small since the 
sensitivities of all selection channels are similar. In the case of dn-type squarks (X'n^ 7^ 0), the 
relative contributions of the decay channels and the resulting limit on A'^j^^ are almost indepen- 
dent of the slepton mass, since gauge decays of Qj^ squarks via charginos are suppressed. 

The branching ratios to the various decay channels depend on the SUSY parameters. Thus, 
for parameter values different from those discussed above, different decay channels are domi- 
nant. For instance, for a zino-like Xi the dominant channels at lower squark masses are those 
with a u and several jets in the final state [10]. Cascade decays of q^ squarks are also possible 
for some parameter configurations via gauge decays involving neutralinos or gluinos. 

6.2.2 Scan of the parameter space 

In order to investigate the dependence of the sensitivity on the MSSM parameters, a scan of 
M2 and n is performed for tan P — 6. Again, sleptons are assumed to be degenerate and their 
mass is set to a fixed value of 90 GeV. Other values for Mf and tan /3 lead to very similar 
results. The parameters M2 and fx are varied in the range 70 GeV < M2 < 350 GeV and 
—300 GcV < fJL < 300 GeV. Parameter sets leading to a scalar LSP or to LSP masses below 
30 GeV are not considered. The latter restriction, as well as the lower value for M2, are mo- 
tivated by the exclusion domains resulting from x searches in I/tp SUSY at LEP [4]. Upper 
bounds on the couplings A'^ ^ and A'^^^^ are obtained for each point in the M2) plane. The 
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Figure 8: Exclusion limits (95 % CL) on A'^^^ for (a) j = 1,2 and (b) j = 3 as a function of 
the squark mass from a scan of the MSSM parameter space as indicated in the figures. The 
two full curves indicate the strongest and the weakest limits on A' in the parameter space in- 
vestigated. Indirect limits from neutrinoless double beta decay experiments (PPOu) and atomic 
parity violation (APV) are also shown. 

results are shown for A'^^-^ in figure 8 and for A'^^^^, in figure 9. For each plot, the two full curves 
indicate the strongest and weakest limits obtained for A' in the parameter space investigated. 
As can be seen from the narrowness of the region that is excluded in only part of the parameter 
space, the limits on both and A'^^^^, are widely independent of the SUSY parameters. For a 
Yukawa coupling of electromagnetic strength, i.e. X'^j^ = y/Ana^ = 0.3 (A'^j^;, = 0.3), m^, cl 
and ii (cIr, sr and b^) squarks with masses below ~ 275 GeV (280 GeV) are excluded at the 
95% CL. For a coupling strength smaller by a factor of 100, masses up to ~ 220 GeV are ruled 
out. 

In figures 8 and 9 the results for the direct production of squarks are compared with indirect 
limits from virtual squark exchange in low energy experiments [29]. The production of u and d 
squarks via a X'm coupling is tightly constrained by the non-observation of neutrinoless double 
beta decay (PPOu) [30]. The best indirect limit on the couplings \[2i and A'^3^ comes from 
atomic parity violation (APV) measurements [29, 31]. The best indirect limit on the couplings 
\[^2 ^iid A'^]^3 results from tests of charged current universality (CCU) [32]. The HERA results 
improve the limits on A' for squarks of the second and third family (i.e. \'i2ii -^'isi; -^'112) -^'113) 
for masses up to ~ 255 GeV. 
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Figure 9: Exclusion limits (95% CL) on A'^^^^ for (a) k = 1,2 and (b) = 3 as a function 
of the squark mass from a scan of the MSSM parameter space. The two full curves indicate 
the strongest and the weakest limits on A'. Indirect limits from neutrinoless double beta decay 
experiments (PPOv) and tests of charged current universality (CCU) are also shown. 

6.3 Limits in the minimal Supergravity model 

In this section the minimal Supergravity (mSUGRA) model [33] is considered, where the num- 
ber of free parameters is reduced by assuming, in addition to the GUT relation between Mi, 
M2 and M3 mentioned previously, a universal mass parameter mo for all scalar fields at the 
GUT scale. By requiring in addition Radiative Electroweak Symmetry Breaking (REWSB) the 
model is completely determined by mo, mi/2, tan/9, the sign of fi and the common trilinear 
coupling at the GUT scale Aq. The modulus of /i is related to the other model parameters. The 
program SUSPECT 2.1 [34] is used to obtain the REWSB solution for and calculate the full 
supersymmetric mass spectrum. 

Assuming a fixed value for the couplings A'^^^ and A'^^^, constraints on the mSUGRA 
parameters can be set, for example on (mo, mi/2), when tan/3, Aq, and the sign of ^ are fixed. 
Aq enters only marginally in the interpretation of physics results at the electroweak scale and 
it is set to zero. The efficiencies for the detection of all gauge decays of squarks involving a 
gaugino lighter than 30 GeV are set to zero since the parameterisation of the efficiencies is not 
valid in this domain. The corresponding parameter space is already excluded by x searches in 
^pSUSYatLEP [4]. 
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Figure 10: Excluded regions (95 % CL) in mSUGRA with A^^ = 0.3 for (a) tan/3 = 2 and (b) 
tan/9 = 6. The region marked "not allowed" corresponds to values of the parameters where 
no REWSB solution is possible or where the LSP is a sfermion. The dashed lines indicate the 
curves of constant squark {ul, ii) mass. The limits from LEP and the TeVatron are given by the 
dotted lines. 

6.3.1 Results for the first and second families 

For fi < 0, the exclusion limits at the 95 % CL obtained for a Yukawa coupling A'^^^ = 0.3 
{j = 1, 2) in the (mo, mi/2) plane are shown by the hatched histograms in figure 10 for the two 
example values (a) tan/? = 2 and (b) tan/3 = 6. The corresponding results for \[^f^ = 0.3 
(A; = 1, 2) are shown in figure 11. The domains marked "not allowed" correspond to parameter 
values where no REWSB solution is possible or where the LSP is a sfermion. 

The constraints on (mo, mi/2) are very similar for both values of tan /3 and both_^p coupling 
types, and A'^^^^, since the mixing of the squark states is very small for j = 1,2 and k = 1,2. 
The excluded regions approximately follow curves of constant squark mass. For X[j^ = 0.3, 
the parameter space defined by Mg < 275 GeV is nearly fully excluded. For A'^^^;, = 0.3, the 
squark mass limit is slightly higher, Mg < 285 GeV, because of the higher squark production 
cross section in e^p collisions for equal couplings. 

The results of the searches for SUSY by the DO experiment [6] at the TeVatron, which 
exploit di-electron events, are also shown in figures 10 and 11. For tan/3 = 2, the HI limits 
are more stringent only for low values of mo, whereas for tan/3 = 6 the domain excluded by 
HI extends considerably beyond the region ruled out by the DO experiment. For tan/3 = 2, 
the parameter space is more strongly constrained by the searches for x's and sleptons at the L3 
experiment [4] at LEP, as shown in figures 10 and 11. This is the only tan/3 value considered 
in [4]. Results for higher values are expected to be similar. The LEP and TeVatron limits are 
independent of the Yukawa coupling. 
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6.3.2 Results on stop and sbottom production 

A non-vanishing coupling would lead to the production of a stop squark. The weak stop 
eigenstates and mix through an angle 9^ to form the mass eigenstates ii = cos O^tL + 
sin 9f iji and ^2 = — sin 9f ii + cos 9f Ir, whose production cross sections scale as A'j^^ cos^ 9f 
and A'j^I^ sin^ 9f, respectively. Thus the lighter state ii does not necessarily have the largest 
production cross section. Similarly, for a non-vanishing X[^^, sbottom production could be 
possible. The weak sbottom states 6^ and mix to form the mass eigenstates 61 = cos 9ihL + 
sm9ibB and 62 = —sm9-^bL + cos^^j6r and the production cross section for 61 (62) scales 
as X'-^i^ sin^ 9^ (X'u^ cos^ 9i). The treatment of stop production is described in the following. 
Sbottom mixing is treated in the same way. 

For the selection channels where the signal is integrated over the whole mass range, the 
fraction of the visible signal in a given selection channel k, is 2(^/^)fc,iC''i/c"tot> where 
{eP)k,i is the total visible branching ratio^ of the selection channel k for the state ii, Ui is the 
production cross section of U and cTtot = (J1+C2 is the total signal cross section. For the channels 
in which the signal is integrated over a mass bin only the contribution of the state U for which 
the sensitivity is maximal, i.e. for which Ci(^fc(/5£)fc,i) is maximal, is taken into account in 
the above summation. The numbers of observed and expected events are then integrated in the 
mass bin corresponding to ii only. 

Using this procedure for both the stop and the sbottom case, exclusion limits are derived for 
= and < for tan P — 2 and tan P — 6. The excluded regions in the (mo, mi/2) plane 
for A'lgi = 0.3 and A'^^g = 0.3 are shown in figures 10 and 11, respectively. The domain below 
the line mi/2 ~ 10 GeV is not considered since it corresponds to cases where the only possible 
LSP decay, into ubd, is kinematically forbidden. 

In the case of stop production for tan /5 = 2, shown in figure 10 (a), the excluded domain is 
slightly larger than that ruled out previously for A'^^ = 0.3 (j = 1, 2) due to the mixing in the 
stop sector which leads to ii masses smaller than the masses of the other squarks. Consequently, 
larger values of mi/2 and mo can be probed. As shown in figure 10 (b) this remains the case 
for tan P — 6 as long as mi/2 is large enough to ensure that the mass of the lightest neutralino 
is above 30 GeV. When the Xi becomes too light, the detection efficiencies for the channels 
involving a Xi (in particular the process xt ~^ Xi) ^re set to zero and the sensitivity is only 
through the eq channel or the decays i bxi followed by a ]/lp decay of the chargino. For 
even smaller mi/2, if the Xi mass is below 30 GeV, only the eq channel contributes. For 
tan P — 2 (tan P — 6), ii masses up to 265 GeV (270 GeV) can be excluded for \[^^ = 0.3. 
These masses are smaller than the maximal sensitivity reached for the same coupling value for 
j = 1,2 because of the cos^ 9f reduction of the ii cross section. For tan/? = 2, the limits 
obtained from this analysis are comparable to the LEP sensitivity in x and slepton searches at 
intermediate values of mo. In the same part of the parameter space, the HI limits for higher 
values of tan P extend considerably beyond the LEP sensitivity which is expected to be similar 
to that for tan /3 = 2. 

*The total visible branching ratio (e/3)fc of a selection channel k is given by (e/3)fe = ^kjf^j, where f3j is 
the branching ratio of the squark decay mode j and j is the corresponding efficiency in the selection channel k. 
The sum runs over all decay modes j considered in the selection channel. 



23 



^ 180 
> 

^160 
"c!l40 
S 120 
100 
80 
60 
40 
20 




mSUGRA, „k = 0.3, tanP=2 



, 1 1 

n < 0, A J = 


1 1 excluded for k=3 


: (a) \ 


\^;^ excluded for k=l,2 : 

HI 




>'''"^\. L3 liinit 

DO limit 

. . 'M; - 






not allowea** 

1,1 





50 100 150 200 250 300 350 

mo (GeV) 



^ 180 
> 



mSUGRA, = 0.3, tanp=6 



H < 0, A „ = 1^ excluded for k=3 
excluded for k=l,2 

HI 




100 150 200 250 300 350 
mo (GeV) 



Figure 11: Excluded regions (95 % CL) in mSUGRA with A'^^ = 0.3 for (a) tan/5 = 2 and (b) 
tan/9 = 6. The region marked "not allowed" corresponds to values of the parameters where 
no REWSB solution is possible or where the LSP is a sfermion. The dashed lines indicate the 
curves of constant squark mass. The limits from LEP and the TeVatron are given by the dotted 
lines. 

In the case of sbottom production (figure 11) the limits are similar to those obtained for 
k = 1,2, since the mixing in the sbottom sector is small at the tan P values considered. Thus, 
the mass difference between bi and 62 is small. The sensitivity follows curves of equal 62 masses 
because the production cross section for this state is much higher than for bi if the mixing angle 
is small. The parameter space leading to 62 rnasses less than 280 GeV is ruled out. 

6.3.3 Dependence of the results on tan /3 

In order to extend the parameter space to larger values of tan /3, a scan of this parameter is 
carried out. The number of free parameters is reduced by setting the masses mo and mi/2 to a 
common value M. The 95 % CL limits on M are shown in figure 12 as a function of tan P for 
^'ijfc — AH squark flavours are considered. For the first two families the exclusion curves 
are rather flat since mixing effects are very small. Assuming equal ^.p couplings, a larger part 
of the parameter space is excluded for d and s production than for u and c production because 
of the higher squark production cross section in e~p collisions. For squarks of the third family, 
mixing effects become important. For tan/3 ^ 10 the increase of the mixing angle 9^ results 
in an improvement of the sbottom limit since it leads to a smaller bi mass, giving a higher 
bi production cross section. The mixing effects are largest in the stop sector, leading to more 
stringent limits on M. For very low values of tan P, the cos^ 9^ reduction of the ii production 
cross section is important. At values of tan /3 ^ 37, the mixing of the two stau (f) states leads to 
decay chains involving light fs which result in final states including r leptons. These channels 
are not searched for explicitly. Thus, in this region of the parameter space, the limit on stop 
production becomes less restrictive. 
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Figure 12: Exclusion limits for mo = mi/2 = M in mSUGRA as a function of tan/?. The 
95 % CL exclusion limits for A'^^ = 0.3 are shown. The areas below the curves are excluded. 
The region marked "not allowed" corresponds to values of the parameters where no REWSB 
solution is possible or where the LSP is a sfermion. 

7 Conclusion 

A search for the production of squarks in e~^p and e~p collisions at HERA has been pre- 
sented. No significant deviation from the SM is observed in any of the final state topologies 
resulting from direct or indirect Rp violating squark decays. Mass dependent limits on the 
couplings are derived within a phenomenological version of the MSSM. The existence of 

M^-type and d^-type squarks of all three generations with masses up to 275 GeV and 280 GeV, 
respectively, is excluded at the 95% CL, for a Yukawa coupling equal to V^tFo^, in a large 
part of the MSSM parameter space. These mass limits extend considerably beyond the reach of 
other collider experiments. For lower squark masses, the results improve the indirect bounds set 
by low-energy experiments. Exclusion limits are also derived in the more restricted mSUGRA 
model, for which the limits obtained are competitive with and complementary to those derived 
at the LEP and TeVatron colliders. 
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